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Abstract 

Chikungunya virus is positive-sense single-stranded RNA virus that causes an 
arthropod-borne chikungunya fever, myalgia and arthralgia. Chikungunya 
virus belongs to the Togaviridae family, and the genus is Alphavirus. Virus- 
host protein interaction plays a vital role in developing vaccines and antiviral 
drugs. We designed the current study to establish the in-silico interaction of 
non-structural protein 2 (nsP2) with proteins of innate immune pathway. The 
nsP2 sequences of various Chikungunya virus genotypes were retrieved from 
National Centre for Biotechnology Institute (NCBI). The homology models of 
proteins were generated through a protein modeling online web server. Protein- 
protein interaction (PPI) between nsP2 and proteins of innate immune pathway 
were docked using High Ambiguity-Driven Docking (HADDOCK) webserver. 
The interactive residues of the bimolecular complexes were analyzed with 
PDBsum-Generate online webserver. Our findings revealed differentially 
affinity of nsP2 of various chikungunya genotypes towards key proteins of 
cellular innate pathway. The nsP2 of Asian genotype demonstrates relatively 
high interaction with interferon-beta promoter stimulator 1 (IPS-1). Similarly, 
nsP2 of various genotypes binds with differential affinity to tumor necrosis 
factor receptor-associated factor 6 (TRAF6) with the highest affinity observed 
for the nsP2 of the West African genotype. Bimolecular complexes of nsP2 and 
host proteins demonstrate the interaction of various domains of nsP2 with 
proteins of the innate immune pathway. Thus, it is sought that the selected 
panel of the proteins might be helpful to treat the viral infection as a therapeutic 
drug target in the future. 
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Introduction 


Chikungunya virus (СНІКУ) is an arthropod-borne 
human pathogen causing chikungunya fever 
associated with arthralgia, myalgia, headache, skin 
rash, and joint swelling [1]. 11953, CHIKV was a first 
time isolated from a patient in Tanzania [2]. CHIKV has 
three genotypes: Asian, West African (WA), and East 
Southern African (ECSA) [3]. СНІКУ belongs to the 
family Togaviridae and the genus alphavirus. 
CHIKV is an enveloped single-stranded positive- 
sense RNA virus having a 5‘methylation cap and 
3*polyadenylation tail. CHIKV genome comprises 
two open reading frames that encode two polyproteins 
and subsequently split into three structural proteins; 
capsid (C) and envelop glycoproteins (E1 and E2) that 
together facilitate viral replication in the host cell. In 
contrast, the four non-structural proteins (nsP1-4) are 
most significant for RNA replication because of their 
enzymatic functions [4,5]. The nsP1 regulates RNA 
capping, nsP2 has helicase and protease activity, nsP3 
is involved in RNA replication, and nsP4 facilitates 
RNA-dependent RNA polymerase functions as shown 
in (Fig. 1) [4,6]. 

Among the СШКУ proteins, nsP2 is responsible for 
the host transcriptional shut-off and facilitates viral 
replication [7]. The nsP2 is five domain long-chain 
protein consisting of 798 amino acid (aa) residues. 
The N-terminal site of nsP2 encodes three domains of 
(1-167 aa residues), while two Rec A domains encode 


535 aa 798 aa 530aa 530 аа 
(-)strand Helicase RNA RdRP 
RNA synthesis Protease synthesis 


261 аа 


Capsid 
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with (168-470 aa residues). Rec A domain is 
responsible for nucleoside triphosphatase and RNA- 
dependent triphosphatase functions [8]. While the C- 
terminal site has two structural domains; the protein- 
like protease (PLP) domain encodes (471-605 aa 
residues), and the methyl transferase-like domain 
(MTL) encodes (606-798 aa residues) respectively. 
The MTL domain is considered non-functional due to 
the lack of significant structural elements. On the 
other hand, the PLP domain has a proteolytic function 
that processes polyproteins and is required for viral 
replication [9]. The CHIKV nsP2 receptor binding 
domain is an ideal site to be the target in the drug 
development pipeline.[10]. 

The current study involves implementing 
bioinformatics tools to predict protein-protein 
interaction of CHIKV nsP2 and human innate 
pathways proteins. The bioinformatics approach relies 
on the protein structure information to determine 
protein interaction with their target proteins [11]. 
Previously, researchers used various bioinformatics 
tools for predicting virus-host interactions for 
obtaining information from protein structure [12,13]. 
Researchers currently using sequence and homology 
models to identify virus-host endogenous and 
exogenous protein interactions [14]. Recently, the 
incredible developments in bioinformatics tools make 
it possible that protein-protein interaction enhances 
the innovation of many potential inhibitors against the 
pathogen [15]. 


ЛЕ Poly A-3’ 


64aa 423аа 61 аа 435aa 
Major envelope (Е) glycoprotein: 
Е1 & E2, Small accessory peptide: 


E3 & 6K 


Fig. 1: Organization of Chikungunya virus whole genome. The figure shows both nonstructural, and structural proteins, 


and 5‘ апаз’ untranslated regions 


Recently, the re-emergence of arboviruses, including 
CHIKV has raised public health concerns because of 
the virus's contagious nature [2]. The arboviral 
diseases have caused large-scale outbreaks in the last 
twenty years, specifically in an endemic country like 
Brazil, which has resulted in high rates of mortality 
and morbidities. Similarly, the emergence of SARS- 


CoV-2 in these arboviral endemic areas was initially 
misdiagnosed because of the coinfection with either 
СНІКУ or dengue [16]. There is no well-developed 
controlling and management strategy existing for 
individuals suffering post-CHIKV consequences for 
many years. This neglect might be because both 
СНІКУ and dengue have been considered a simple 
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and short-lived infection [17]. Therefore, it is more 
desirable to establish controlling and therapeutic 
strategies against CHIKV. 


Materials and Methods 


Protein sequences retrieval and modeling 


We retrieved nsP2 sequences of CHIKV genotypes; 
Asian, WA, ECSA, and human innate immune 
pathway proteins from the NCBI. Through SWISS- 
MODEL webserver 
(https://swissmodel.expasy.org/interactive), we 
modeled the Three-dimensional (30) structure of 
CHIKV nsP2 and proteins of the human innate 
pathway such as; retinoic acid-inducible gene I (RIG- 
I), interferon-beta promoter stimulator 1 (IPS-1), 
Tumor necrosis factor receptor-associated factor 3 
(TRAF3), Tumor necrosis factor receptor-associated 
factor 6 (TRAF6), and ТІК domain-containing 
adaptor inducing interferon-B (TRIF).SWISS- 
MODEL first aligns the amino acid sequence of the 
target protein and uses this sequence as a query against 
the template library. Further select the best template 
on the base of sequence identity, estimated quaternary 
structure, and global model quality estimate then 
generates а 3D structure of the protein. 
The stereochemical stability of the modeled proteins 
was analyzed through the ERRAT program. The 
higher ERRAT scores represent good qualities of 
protein. Generally, a score above 50 was considered a 
good quality model [18]. The ERRAT scores of both 
virus and host were shown in Table S1. 


Docking of nsP2 and Proteins of Innate Immune 
Pathway 


The energy of proteins was minimized with molecular 
operating environment (MOE) software, using the 
Amber 10 force field by adding missing hydrogen and 
partial charge to proteins [11]. To gain insight into 
protein interaction, we docked nsP2 and innate 
immune proteins through an online web server, High 
Ambiguity Driven protein-protein 
DOCKing (HADDOCK), and = pyDockWEB [19, 
20,21]. 


Analysis of biomolecular complex and interacting 
residues 


We analyzed the bimolecular complex using UCSF 
Chimera (version 1.13.1) [22]. The PDBsum- 
Generate(http://www.ebi.ac.uk/thornton- 
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srv/databases/pdmsum/Generate.html) online tool 
was used and determined various kinds of interactions 
between the interactive residues of virus-host 
proteins. 


Results 


HADDOCK webserver provides the top ten 
biomolecular complexes. We selected the best 
complex based on the HADDOCK score, 
electrostatic. and Van der Waals energies. 
HADDOCK score is the average of all given energies 
[23]. The statistical scores of the HADDOCK 
complexes were shown in Table 1. Best complexes 
have the lowest electrostatic, van der Waals energies, 
and HADDOCK score. The stability of two proteins is 
based on the number of active residues involved in the 
complex formation and shows the interface between 
the two molecules. Large interfaces between the 
molecules show high binding energies [23]. We also 
used the pyDockWEB server to validate the 
HADDOCK complexes. pyDockWEB provides both 
protein-protein/ligand interaction to the user. 
pyDockWEB generated several complexes [20]. The 
best complex was selected through their given 
energies such as electrostatic, desolvation Vander 
wall, and total energy, as shown in Table S2. 

The sequence alignment data shows that CHIKV nsP2 
was highly conserved within genotypes. There were 
very few mutations among all genotypes as shown in 
(Fig. 2).The PDBsum-Generate results showed 
various kinds of interaction involved in the stability of 
protein-protein interaction. CHIKV nsP2 showed 
various numbers of hydrogen bonds, and salt 
bridge interaction for either genotype-specific/host, 
as we described in Table 53-97 and Fig. 52-57. 
HADDOCK complexes have various levels of 
interaction affinity between nsP2 of different CHIKV 
genotypes and various cellular proteins of the innate 
pathway. RIG-I senses single-stranded RNA viruses 
to induce antiviral responses. In the current study, 
RIG-1 was included as a positive control to determine 
protein-protein interaction among CHIKV nsP2 and 
various innate immune proteins. The HADDOCK 
score of RIG-1 and CHIKV genotypes, Asian, WA, 
and ECSA, were -103.7 + 8.4, -105.7 +7.2, and -108.5 
+ 2.3. HADDOCK score of the biomolecular 
complexes of nsP2 of СШКУ all genotypes and the 
human innate immune proteins has comparable 
binding affinity with the RIG-1 were described in 
Table 1 and complexes were shown in Fig. S1A. The 
pyDock complexes score were provided in 
supplementary material as described in Table 52. 
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Table 1: Analysis of HADDOCK complexes of nsP2 of CHIKV genotypes and human innate immune proteins 


СНІКУ genotypes -host HADDOCK Score 


Van der Waals energy 


Electrostatic energy 


Desolvation energy Restraints violation 


proteins (a.u) (Keal/mol) (Keal/mol) (Keal/mol) energy (Kcal/mol) Z-Score 

Asian- nsP2- RIG-1 -103.7 + 8.4 -86.4 = 7.0 -260.2 + 44.8 -26.4 = 5.2 611.4 = 69.95 -1.5 
WA- nsP2- RIG-1 -105.7 = 7.2 -91.2 = 4.6 -364.4 + 37.0 -9.1 + 2.9 675.4 + 113.62 -1.3 
ECSA-nsP2- RIG-1 -108.5 + 2.3 -85.6 =4.6 -392.8 + 35.4 -8.1 = 7.4 637.2 + 59.58 -1.6 
Asian- nsP2- IPS-1 -120.144.4 -89.1 = 4.3 -332.6 + 49.4 -23.6 =2.9 591.2 + 2.86 -1.7 
WA- nsP2- IPS-1 -101.0+6.1 -75.3 = 15.7 -370.6 = 88.7 -95.8 = 13.8 1441.5 = 103.17 -2.1 
ECSA- nsP2- IPS-1 -114.347.9 -73.5 = 8.4 -345.2 + 16.0 -33.6 + 1.8 618.6 + 111.82 -1.7 
Asian- nsP2- TRAF3 -102.9 + 7.4 -89.5 + 2.8 -352.5 + 56.1 -36.1 + 7.8 931.8 + 109.17 -1.5 
WA- nsP2- TRAF3 -110.1 + 16.2 -97.3 + 1.7 -377.4 + 22.9 -9.6 + 8.0 722.7 + 113.14 -1.9 
ECSA- nsP2- TRAF3 -113.3 + 4.5 -89.1 = 6.2 -338.9 = 57.9 -38.4 = 3.9 819.8 + 167.76 -1.9 
Asian- nsP2- TRAF6 -102.1 = 4.3 -72.9 + 6.4 -419.9 + 81.7 -0.6 + 9.5 553.7 + 44.49 -1.2 
WA- nsP2- TRAF6 -93.9 + 8.1 -71.6 15.0 -302.3 + 55.6 -24.1 = 8.2 622.1 + 71.88 -1.6 
ECSA- nsP2- TRAF6 -107.4 =2.1 -71.3 = 4.9 -437.9 + 46.7 -0.6 = 8.8 521.5 + 27.50 -1.4 
Asian- nsP2- TRIF -117.6+5.6 -89.5 = 2.8 -218.9 + 45.7 -31.0 = 8.8 467.7 + 124.48 -1.5 
WA- nsP2- TRIF 112.6 = 6.9 -91.5 = 4.5 -212.1 + 42.3 -35.9 = 7.8 572.6 + 47.01 -1.9 
ECSA- nsP2- TRIF -113.5 = 6.6 -95.0 = 2.4 -322.2 + 61.2 -19.0 = 6.6 649.4 + 18.45 -1.3 


Note: arbitrary unit (a.u), Kilocalories per mole (kcal/mole), HADDOCK is the sum of all energies involved in protein-protein interaction, Z- 


score is the standard deviation of the energies 


Among the HADDOCK complexes, the nsP2 CHIKV 
Asian genotype with IPS-1 showed comparatively 
high binding affinity while nsP2 of WA genotype with 
TRAF6 showed low binding affinity -93.9 + 8.1. Both 
CHIKV nsP2 and host innate immune pathway 
proteins were labeled with different colors for better 
understanding. In the HADDOCK complex, host 
cellular proteins were presented by green and CHIKV 
nsP2 in red and blue colors. The red color represents 
the PLP domain, while the blue represents the MLT 
domain of nsP2 of CHIKV (Fig. S1 A-E). 

In the biomolecular complexes chain, A represents 
CHIKV nsP2 and chain B represents host immune 
protein. The human host innate immune proteins 
showed various interactions among CHIKV 
genotypes, such as RIG-I and nsP2 of Asian genotypes 
with one disulfide bond, 15 hydrogen bonds, and 189 
non-covalent bonds. Similarly, the WA genotype 
shows one disulfide, 16 hydrogens, and 213 non- 
covalent bonds, whereas the ECSA genotype shows 
one disulfide, 21 hydrogen, and 230 nonbonding 
interactions. The biomolecular complex of IPS-1 and 
СНІКУ genotypes; Asian, WA, and ECSA showed 4, 
2, 2 disulfide bonds, 16, 14, 16 hydrogen bonds, and 
149, 184, 180 non-bonding interactions, while 
CHIKV nsP2 -TRAF3 complex showed 4, 3, 2 
disulfide bonds, 19, 14, 1 hydrogen bonds and 211, 
169, 183 non-bonding interactions respectively. The 
TRAF6 and TRIF interaction with CHIKV nsP2 
showed various interactions between Asian, WA, and 
ECSA genotypes. TRAF6 showed 7, 4, 4 disulfide 
bonds, 17, 12, 14 hydrogen bonds and 160, 146, 186 
non-bonding interactions, whereas nsP2-TRIF 
complex showed 0, 1, 1 disulfide bonds, 19, 16, 21 
hydrogen bonds, and 179, 213, 230 non-bonding 


interactions. The disulfide, hydrogen, and non- 
bonding interactions showed in red, blue, and orange 
colors (Fig. S2 A-E). The interactive residue involved 
in the disulfide and hydrogen bond formation was 
described in Table S3-S7 and Fig. S3-S7. 


Discussion 


The nsP2 is the largest protein among the nsPs of 
CHIKV. Based on the helicase and protease function, 
the nsP2 of СШКУ is targeted for molecular docking 
to achieve а significant result. Protein-Protein 
interaction at specific amino acid residues is 
considered essential in molecular docking. The 
bioinformatics approaches were considered a fast and 
reliable method to predict the potential interaction 
between virus-host cellular proteins, which played a 
significant role in developing antiviral drugs. The 
Three-dimensional (3D) structure of protein provides 
an opportunity to dock various immune cellular 
proteins against viral proteins [11,12]. Previously 
researchers conducted a study on docking-based 
interaction of proteins, an antimicrobial compound, 
and phytochemicals, and in vitro and in vivo 
experiments are still going on. Protein-protein 
docking has become an essential tool for antiviral drug 
discovery. Hence, the current study concentrated on 
the bioinformatics approach toward evaluating the 
antiviral effect of these six innate immune proteins on 
СШКУ nsP2. The СШКУ nsP2 was used as a target 
protein for these immune proteins to identify potential 
antiviral drugs for CHIKV infections [24]. 

The innate immune system provides the first line of 
defense in response to microbial infection. Pathogen- 
associated molecular patterns (PAMPs) are the 


165 


Asain Е | НЯ | |. тины McEREY@EsP В р ah 
co ЕЕ ЕЕ НИ: 
WA cum ITPRGA Н Е 


болели, СГ ТЕТРЕСА | СИЕ VGEYLVLSPQ Ерт 


СМА 


Biomedical Letters 2022; 8(2):162-169 
BAB инсвхс Ши AMBARGCREE МР 
СМА 


80 
1 
Е] 
ЕЕ СТ ЕИ Е В 
HsG Г Мб с РБС АШ: 


Conservation. 


m ПЦ 


a 


НН 
a B 


RRERAMAmAG 
A BEQSESESAT RREAAMABAC РАШМТОЕ 
сашкан. DFQSLSESAT М 


УМ RKLHHIAMHG PALNTDE 


Аза" МНАВАМЕ с РА [<] 


90 
ВК ШЕРАСРИК Bcuece колк 
EGSA + EAHECER ШЕРА eksa 


Cons: ervaton. 


|: рор 

PERMA В сМРС 
MBHEBAMECER ШЕРАСРИКТА МИСИЙСМРСЗ СИЗАШИК 

aowa, YHEFAYEGLK IRPACPYKIA VIGVFGVPGS GKSAI IK 


Conservation 


Е 
E 


1 


Esi BEURAER 


PP 160 


іі амаз CTHSGRAGRY AVEAYDGRVL пов 
Г Ня Е 160 
60 


M TT Ш] 


ны 160 
H Е рир в ККВЕлАСШ ЕЕ 
ESY ERO 


CRREBAAGCH ЕАН 
ПА О О А TEE 


700 “ne 20 
a = He 


GK Г: p И т тава at 


NLV TRQDLVTSGK КЕМСОЕ | ТТО VMRQRGLE!S ARTMDSILLLN 


240 
240 
240 


Азап ССМИРМОМШМ МОБАВАСН$С НА 


“Wa CONRPMOMEN HBEABACHSc THUAUHAGUR 
ВР И MBEABACH ABMABYR P 


МА G 


Consensus. GCNRPVDVLY VDEAFACHSG TLLALIALVR PRQKVVL 


PRORMMECGE На a 


АЕ PROCCHENMM © 


| | Ш ШИ L 
н 


ae 
uuananEù OMMHKsiis BRCTECETAR 320 
CGD PKQCGFFNMM QMKVNYNHNI 


П ПЦ 


Cons: eration 


RUNENHNH CTOMMAKSMS RACTEPRT All 520 
ШО 


BEL 320 
BROT eC eat 


Ш 


| 
F 


SAE 


до Е ИНМЕСК ee HBr Ice rk 
CSA G 


Ес 
Consensus УЗЗЕНУЕСКМ TU е тте 


CTQVYHKSIS 
380 


390 400 


1 
той Row В ЕСЕ ва Cl Ha HH 400 
E acum ШО ЕСЕН cETRKCE MANR EN 400 
[ег te Te] ERGHEMMTAA ASQGGHETRKGY МАМНОКЫМЕМ 400 


ТСЕ ЕНЕ YRGYEVMTAA ASQGLTRKGV YAVRQKVNEN 


EWA assume GH 


Cons: ervaton 


ни 


ПШ ШЕТШ 


90 


Азап PEMASTSEAW RTECH BUWRTE GBP WERTE 
ЕСЗА PEMASTSEHM i Еск BEWK i cBP WiK 
wa PBRASTSEHM BTEcK BMWKTEscBe ШКТ 


Consenets: PLYASTSEHV RE LVWKTLSGDP WIKTLQN 


надр PT Ш 


so 


Conservation 


440 2160 


1 Ka 
PPK СМИКА HEHA H ЕВ. wa 480 
PPK GNEKA Ба BRENA а WA 480 
PPK СМИКА НЫ ЕМЕНА А 480 


PPK СМЕКАТ | KEW ит PTET TE 


‘ii al cca нЕ 


an i 


50 


Аѕап ИЗШМРИШЕ! A о MGABKEBK ARSPERA 0-8 ВЕ con РЕЙ 07] НАЗЕН ЕСЕ 560 
Econ | НИ | AH | all НЯ i Tava cHESKPEY HHI HE | RPGGKMEGEH 560 
ЕРШВЕ' Г CHESKPEY SHAMA PGCKMEGE 550 


белам, CRNE $ rkt nonas QI ЕЕ КОРЕЯ 


Conservation 


ine ICTRMYGVDL DSGLFSKPLV SVYYADNHWD NRPGGKMFGF 


| Ц 


Азат 
ЕСЗА 
WA 


| Е Ш 
В : Е Я пава 
АА iN] 
P 


ШЕЕ КУРЕ E |. ЕВЕ 
NPEAASILER KYPFTKGKWN TE TE 


WER КИРЕТ RHEDEN? 


ИРЕ 


CONSENSUS. 


600 


1 
H nai НЫЕ свв 640 
НЕ $ RMEWENNK MNGHHMEBMS 510 
N ANRREP ABH ENK BNCHAME 640 
TTN EEN SLVAEHRPVK ТТ Т 


снаа, ШП 11 AE 0 Ш = Ш 
0% 
seo Ean 00 zza 
Asan GHNNAHPTER H APEcHE САШИ BABES cHPATECREO HUUNNMETPE винна р Haul EREERE SC 720 
ЕСЗА APBCHR GABNTENBER сШРАТИСВЫО EMWENEHTPR ВИННИ PGG 720 
wa GHNEMEPTER MTWHAPECHR ЅАВИТИМИВИ CHPATECHNE ШИВИНИНТРИ BEHHEG а HANK EREEKPSS 720 


GYNLALPTKR VTWVAPLGVR 


ее 


Consens ge 
Я 


n || 


Conservation 


HRPP 
HRPP 


Herel 


LKPPCVT 


| 


ia 
| ABA GR EcREER 
HEA OR BcRAER 


Вис ESRA 
ВЕШВАМСИл ВЕ EcRKER 


8, SLLIRAYGYA ОКТЗЕКУ | СУ LGRKFRSSRA 


Asain 


ganean 
ns 


EEI 


азн 


a 
Е au |=. 
BTTHY СОА 
Еттн@ сал 


ТЕО 


7 
ED 
EB 
ЕО 


ВАС 
ВАС 
ВАС 
SNT EMFFLFSNFD NGRRNFTTHV MNNQLNAAFV СОАТВАСС 


798 
798 
798 


GRR 
GRR 
GRR 


ААВ 
ААЙМ 
ААВ 


| Ц 


Fig 2: Consensus sequence of CHIKV nsP2 of Asian, WA, and ECSA genotypes. The consensus sequence has been 
shown by one-letter amino acid codes at the bottom above the bars. The bars show the degree of conserved amino acid 


at their respective position in terms of percentages. 


the microbial substance that pathogen recognition 
receptors (PRRs) recognize on host cell surface 
receptors and activate the innate immune response. 
Thus, the innate immune response prevents or 
eliminates microbial infection. During viral infection, 
host PRRs recognize the viral substance through 
PAMPs and activate the interferon response to inhibit 
viral infection. Toll-like receptors and RIG-like 


receptors are the two main types of PRRs; both 
recognize microbial infection and trigger the immune 
response [25,26]. 

RIG-I detects viral RNA in the cytoplasm, The 
conformational change in RIG-I or МРА-5 is required 
for the interaction between the tandem CARDs of the 
helicases with the CARD of IPS-1, which is 
incorporated into the mitochondrial outer membrane 
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by a transmembrane domain in its C-terminal site. 
Thus, this interaction induces interferon-sensitive 
genes and inflammatory inflammation. Type 1 
interferon (IFNs) plays a significant role in the 
inhibition and clearance of viral infection [27]. In vitro 
study reported that RIG-I could detect synthesized 
RNA through 5’ triphosphate of single-stranded RNA 
[ 28, 29]. In the current study, RIG-I was used as a 
positive control, RIG- I had strong bonding interaction 
with the PLP domain of CHIKV nsP2. The 
HADDOCK scores were relatively similar among the 
CHIKV genotypes; Asian, WA, and ECSA were 
followed (-103.7 + 8.4, -105.7 + 7.2, and -108.5 + 2.3). 
The binding energies of the HADDOCK complex 
were shown in Table 1, and bonding interactions were 
shown in Table S3 and Fig. S3. 

The conformational change in RIG-I or MDA-5 is 
required for the interaction between the tandem 
CARDs of the helicases with the CARD of IPS-1, 
which is incorporated into the mitochondrial outer 
membrane by a transmembrane domain in its C- 
terminal site [22]. While the interaction of these 
CARD-CARD domains, results in the propagation of 
downstream signaling leading to antiviral and 
inflammatory responses [30]. Previously, it has been 
reported that IPS-1 playing important role in antiviral 
immunity [31]. The bimolecular complex of CHIKV 
nsP2 and IPS-1 showed a higher HADDOCK score 
among the targeted host proteins. The HADDOCK of 
all three genotypes; Asian, WA, and ECSA were (- 
120.1 + 4.4, -101.0 + 6.1, and -114.3 + 7.9). The 
binding energies of the HADDOCK complex were 
shown in Table 1, and bonding interactions were 
shown in Table S4 and Fig. S4. 

TRAF3 plays an essential role in the promotion of 
signaling pathways by using interferon regulatory 
factors to produce IFNs which plays role in viral 
clearance [32]. TRAF3 regulates the stimulation of 
TLR-mediated mitogen-activated proteins like TLRs, 
Janus kinase, signal transducer, and activator of 
transcription pathways. These pathways are activated 
in response to viral infection; thus, the cell achieves 
the antiviral state. TRAF3 signal is essential for the 
activation of canonical and non-canonical nuclear 
factor-kappa B (NF-«B) pathways to activate antiviral 
immune response [33]. TRAF6 is one of the recently 
discovered mammalian proteins of the TRAFs family. 
TRAF6 regulates multiple signaling pathways, such 
аз Proto-oncogene tyrosine-protein kinase (Src) 
family tyrosine kinase, an inhibitor of nuclear factor- 
кВ, and mitogen-activated protein kinase. Myeloid 
lineage hematopoietic cells express receptors on the 
cell surface like receptor activators of nuclear factor 
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kappa beta, Cluster of differentiation 40, TLRs, and 
TRAF6, using this for the signaling to stimulate the 
innate and adaptive immunity towards pathogenic 
infection [34,35]. Both HADDOCK score and binding 
energies of TRAF3 and TRAF6 complexes were 
shown in Table 1 and bonding interactions were 
shown in Table 55-56 and Fig. S5-S6. 

TRIF is a TIR-containing protein responsible for the 
activation of the TLR4-dependent pathway and 
Myeloid differentiation primary response 88 
independent pathway, which activates МЕ-кВ and 
IFN-B production. Interaction of TRIF on the C- 
terminal site with receptor-interacting protein 1 is 
essential for the activation of МЕ-кВ. The TRIF 
signaling pathways regulates several transcription 
factors such as interferon regulatory factor 3, activator 
protein 1, and NF-«B, which leads to type 1 INF and 
cytokines production [36]. Both TRIF-dependent and 
independent pathways play a significant role in TLR 
signaling responses and in controlling viral replication 
and pathogenicity. We found a strong bonding affinity 
among the TRIF and nsP2 of CHIKV genotypes. The 
HADDOCK score and bonding interactions were 
shown in Table 57 and Fig. 57. 


Conclusion 


In conclusion, the binding affinity of innate immune 
proteins seems the most potent inhibitor against 
СНІКУ nsP2. Both PLP and MTL domains of CHIKV 
nsP2 were involved in the interaction. Residues 
involved in interaction may be the potential target for 
the drugs and peptide inhibitors. The in vivo approach 
of the selected protein and interacting residues may be 
useful for designing pan-chikungunya antivirals in the 
future. 
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